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In this paper, we propose a novel GL double layer EM cloak in the broad frequency band. In 
short, we call it as GLWF double EM cloak. The GLWF double layer cloak (GLWF cloak) consists 
of two sphere annular layers, two type cloak materials are proposed and installed in its each layer, 
respectively. The outer layer of the GL cloak has the invisible function in broad frequency band, 
while its inner layer has the fully absorption and rapid delay function. The outer layer cloaks 
the Local concealment from the Global exterior EM field; The inner layer cloaks the Global free 
space region from the Local field excited inside the concealment. The GLWF double layer cloak 
overcomes the following difficulties of the single layer PS cloak. (1) There exists no EM wavefield 
can be excited inside the concealment of the PS cloak, the concealment of the PS cloak is blind. 
Our GL double layer cloak recovered that the EM wave field can be excited inside the concealment 
of the GL double layer cloak. (2) There is exceeding light speed physical violation in PS cloak, its 
invisibility only in very narrow frequency band. Our GLWF double layer cloak corrects the violation. 
(3) The reciprocal law is satisfied in our GL double cloak media. However, the PS cloak damaged 
the reciprocal law.. The simulations and comparisons of the EM wave field propagation through the 
GLWF double cloak, GL double cloak and GL double layer with PS outer layer are presented to 
show the advantages of the GLWF double layer EM cloak The GL double layer cloaks are proposed 
by our GL EM modeling and inversion. The 3D GL EM modeling simulations for the double layer 
cloak are presented. The copyright and patent of the GLWF double layer cloak materials and GL 
EM modeling and inversion in this paper are reserved by authors in GL Geophysical Laboratory. 

PACS numbers: 13.40.-f, 41.20.-q, 41.20.jb,42.25.Bs 



I. INTRODUCTION 



We have proved that there exists no Maxwell wavefield 
can be excited by sources inside the single layer cloaked 
concealment with normal materials [1]. For covering this 
difficulty, we proposed a GL double layer cloak [2]. Its 
outer layer has invisibility from the exterior light and 
EM wave field, never disturb the exterior field, and cloak 
the Local inner layer and concealment from the Global 
exterior wavefield. Its inner layer cloak absorbs the in- 
ternal wave field, such that the internal wave field, which 
is excited from the sources inside the concealment, can 
not propagate outside of the inner layer. Our GL double 
layer cloak is proposed by the GL EM modeling and in- 
version [3] . The GL double layer cloak overcomes the fol- 
lowing difficulties of the single layer PS cloak by Pendry 
[4]. (1) The PS cloak damaged the EM environment of 
its concealment, such that there exists no EM wavefield 
can be excited inside the concealment of the PS cloak, 
the concealment of the PS cloak is blind. Our GL double 
layer cloak recovered the normal EM environment in its 
concealment, such that the EM wave field can be excited 
inside the concealment of the GL double layer cloak. (2) 
The two sources reciprocal law is very important princi- 
ple in the electromagnetic theory and application. The 
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reciprocal law is satisfied in our GL double cloak me- 
dia. However, the reciprocal law is damaged by the PS 
cloak. The PS cloak is strong dispersive cloak material, 
in which the e r ee = Rl(r - R2) 2 / 'r 2 / '(R3 - R2) 2 , when 
r — s- i?2 £ r £e ~ (r — R 2 ) 2 0. such that the EM wave 
velocity exceed the light speed in some part of the PS 
cloak that is violated to the physical principle. The PS 
cloak has invisibility only in very narrow frequency band. 
In the GL outer layer cloak [2] [5] [6] , the weak dispersion 
is obtained, when r — > R2 e r ee ~ (r — R2) 1 ' 5 — > 0. In this 
paper, we obtain very significant progress to overcome 
the difficulty, very pefect weak dispersive and degenera- 
tive rate, e r ee s» l/log(r — R 2 ) — > 0, is obtained. Chen 
ct al proposed an analytical method for analysis of the 
PS cloak [7]. 

Finding and exploration is inverse problem; Hiding and 
cloaking is other inverse problem. They have close rela- 
tionship. Based on the 3D GL EM modeling simula- 
tions [3] and GL Metre Carlo inversion [5], We propose a 
novel GL double layer cloak in the broad frequency band. 

The President Professor Yuesheng Li in Sun Yat-Sen 
University very concerns and encourages our research 
works on the GL EM modeling method and the GL dou- 
ble cloak. Our paper is for celebrating his great scientific 
computational carrier in 60 years and his 80 birthday. 

The description order of this paper is as follows. We 
have introduced single layer cloak and our GL double 
layer cloak and main content and research cue path in 
Section 1. In Section 2, we propose a novel GL double 
layer cloak in the broad frequency band The EM wave 
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propagation simulations in our wide frequency band GL 
double layer cloak are presented in Section 3. The im- 
portant reciprocal law in the cloak media is presented in 
Section 4. In Section 5, we will conclude our paper. 



II. THE GL DOUBLE LAYER CLOAK 
MATERIALS IN BROAD FREQUENCY BAND 

For overcoming the weakness of the single layer cloak, 
we proposed the GL double layer cloak in broad fre- 
quency band, in short it is called GLWF, which consists 
of the inner layer cloak and outer layer cloak 



A. GLWF Inner layer Cloak Anisotropic Material 

On the inner sphere annular layer domain, Ugli — 
{r : Ri < r < R 2 } , by the GL EM modeling and inver- 
sion [3] [5], we propose an anisotropic metamaterial as 
follows, 



[ D \gli = dia 9 > 
£i — diag [e r ,i, £0,i, £<t>,i] £b, 
fit = diag \p, r ,i, M0,»> H<^,%[ Mb, 



(1) 



R\-R\ Ri 



Ri 



The SIgli is called as GL inner layer cloak, the materials, 
[D] GLI — diag[si,fii\ in (1), are the anisotropic GL in- 
ner layer cloak metamaterials. Where the subscript GLI 
means the GL inner layer, the symbol diag denotes the 
diagonal matrix, [D]gli is 6 x 6 diagonal matrix, is 
3x3 dielectric diagonal matrix in the inner layer, the 
subscript i denotes the inner layer, \ii is 3 x 3 magnetic 
permeability diagonal matrix in the inner layer, e r ^ is the 
relative dielectric cloak metamaterial which is formulated 
by the fourth sub equation in (1), the subscript index r, i 
denotes the dielectric is in r direction and in the inner 
layer, eg y i is the relative dielectric cloak metamaterial in 
direction and in the inner layer which is formulated 
by the fifth sub equation in (1), The Mr,i is the relative 
permeability cloak metamaterial which is formulated by 
the fourth sub equation in (1), Eb is the basic dielec- 
tric in free space, e b = 0.88541878176 . . . x 1CT 11 F/m, 
Mb is the basic magnetic permeability in free space, 
Mb = 1.25663706 x lQ- & mkg S - 2 A- 2 , other symbols in 
(1) have similar explanation. 

The EM wavefield of sources located inside inner layer 
or the cloaked concealment is completely absorbed by 
the inner layer and never reaches the outside boundary 
of the inner layer, also, the EM wavefield excited in con- 
cealment is not disrupted by the cloak. The inner layer 
metamaterial, in equation (1), cloaks outer space from 
the local field excited in the inner layer and concealment, 



which can also be useful for making a complete absorp- 
tion boundary condition to truncate infinite domain in 
numerical simulation. 



B. GLWF Outer layer Cloak Anisotropic Material 

We proposal a new novel GLWF outer layer cloak in 
this section. Let the outer sphere annular layer domain 
Qglo — { r '■ R2 < t < R3} be the GL outer layer cloak 
with the following anisotropic GLWF outer layer cloak 
metamaterials, 

[D\gLWFO = diag [Sglwfo, figluifo] j 
Sgiwfo = diag [e r , ,glwfo-> E-9,,glwf 07 £<f>„glwfo\ ^bi 
f^glwfo — di(Xg [flr,glwfoi f^9,glwfoi i^<fi,glwfo] f^bi 
£r ,glwfo — ^r.glwfo — 2 (r — R 2 ) 

log (eV^iry^ _ R 2) j ( r 2 _ Rl)/R^/R 2 / r (2) 

£9 : glw fo — \^9 ,glw fo — £<fi,glwfo — f-^<f).glwfo 

«i/(r-ik)/ 

log 1 ' 5 (eV^ry/ifiZ - Rj) / (r 2 - R 2 2 )/R^j /r 
The GL outer layer clok [2] is presented as follows 

\ D ]gLO = dia 9 [£glo,p.glo] , 

Sgio — diag [e r , g io, £e, g io, £<t>,glo] £ b, 

figlo = diag [fJ, r ,glo, V>0,glo, fJ-4>,glo] Mb; 



R 3 r 2 — R 2 \/r 2 — R\ 
^r.glo — f-^r.glo — r 2 2 ~~' 



sj R\-R\ ' 

£@.glo — 1^9 .glo — &cf).glo — f^(p : glo 
_ _gs r 

We chose the PS cloak [4] as outer layer cloak 
\ D \pso = dia 9l 

Sp SO — didg [£r,pso5 £9,psot £</>,pso] £bi 
f-Lpso — didg \Hr.psoi f-^9,psoi f^<p,pso\ i-^bi 

Mr, pso - R 3 (r - R 2 ) 2 /r 2 /(R 3 - R 2 ) 



(3) 



(4) 



-r.pso 



S0,pso — H0,pso — £<f>,pso — H4>,pso — R3K.R3 — R2) 



Where the subscript GLWFO means the GLWF 
outer layer, [D]qlwfo is 6 x 6 diagonal matrix, e g i w f is 
3x3 dielectric diagonal matrix in the outer layer, the sub- 
script o denotes the outer layer, fbgiwfo is 3 x 3 magnetic 
permeability diagonal matrix in the outer layer, £ r , g i w f 
is the relative dielectric cloak metamaterial which is for- 
mulated by the fourth sub equation in (2), the subscript 
index r,glwfo denotes the dielectric is in r direction 
and in the outer layer, e^^/o is the relative dielectric 
cloak metamaterial in <j> direction and in the outer layer 
which is formulated by the fifth sub equation in (2), The 
Vr,glwfo is the relative permeability cloak in outer layer 
and in r direction which is formulated by the fourth sub 
equation in (2), Similar explanation for GL outer layer 
cloak GLO and PS cloak as outer layer PSO, The outer 
layer cloak provides invisibility, does not disturb exterior 
EM wave field, and cloaks the local concealment from the 
global exterior EM wavefield. 
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C. GLWF Double Layer Cloak 

The GL inner cloak Ugli domain and GL outer cloak 
^GLwfo domain are bordering on the sphere annular 
surface r = R 2 . We assemble the Hqli as the inner 
sphere annular domain and Qqlwfo as the outer sphere 
annular domain and make them coupling on their inter- 
face boundary annular surface r = R 2 as follows, 



^GLWF = ^GLI U ^GLWFO 

= {r : Ri < r < R 2 } [j {r : R 2 < r < R 3 } 
= {r:Ri<r<R 3 }, 



(5) 



and situate the double layer anisotropic dielectric and 
magnetic permeability susceptibility tensors [D] GLWF on 
the flcLWF as follows, 



[D] 



GLWF 



[D} GLI ,r en GLI 
[D\glwfo ' r e SIglwfo- 



(6) 



The GL cloak material [D] GLI = diag [si, m] in (1) 
on the inner layer domain £Igli and GL outer layer 
cloak material [D] GLWFO = diag [e g i wfo , p. g i w f ] in (2) 
on outer layer domain SIglwfo are assembled into the 
GL anisotropic doubled layer cloak material on the do- 
main Qglwf- The domain Qglwf with the metamate- 
rial [D] GLWF in (6) is called as the GL double layer cloak 
in broad frequency band, say GLWF double cloak. Sim- 
ilar explanations are for GL double layer cloak coupled 
by GLI in (1) and GLO in (3), and GLI - PSO double 
layer cloak coupled by GLI in (1) and PSO in (4). 

The GL cloak means that the outer layer has invisibil- 
ity, never disturb exterior wavefield and cloaks the Local 
concealment from the exterior Global field, while the in- 
ner layer cloaks outer Global space from interior wave- 
field excited in the Local concealment. The properties 
of the double layer material are studied by GL method 
simulations and analysis in next sections. 



III. THE GL EM SIMULATIONS AND 
COMPARISON OF THE GLWF DOUBLE 
CLOAKS AND PENTRY CLOAK 

A. The Simulation Model of The GLIWFO, GLIO, 
GLI-PSO Double Layer Cloak 

The simulation model: the 3D domain is 
[—0.5m, 0.5m] x [—0.5m, 0.5m] x [—0.5m, 0.5m], the 
mesh number is 201 x 201 x 201, the mesh size is 0.005m. 
The electric current point source is defined as 



6(r - r s )S(t)e, 



(7) 



where the r s denotes the location of the point source, 
the unit vector e is the polarization direction, the time 
step dt = 0.3333 x 10~ 10 second, the frequency band 
is from 0.05 GHz to 15 GHz, the largest frequency 
/ = \hGHz, the shortest wave length is 0.02m. The 



inner layer CIgli is denoted by (1). The EM GL double 
layer cloak Q, G lwf = &gli U ^glwfo in (2), VL G lio = 
^gli U ^glo in (3), and VL G li-pso = ^gli U ^pso in 
(4) are consist of the double spherical annular with the 
center in the origin and interir radius R\ — 0.21m, med- 
dle radius R 2 = 0.31m. and exterior radius R3 = 0.45m. 
The cloak is divided into 90 x 180 x 90 cells. The spherical 
coordinate is used in the sphere r < JS3, the Cartesian 
rectangular coordinate is used in outside Qgl to mesh 
the domain. 



B. Point Source Si In The Concealment And 
Other Point Source S2 In The Free Space 



The two point sources S\ and are used to excite 
the EM wave propagation through the GLWF, GLO, and 
GLI-PSO double layer cloaks. The first point current 
source S\ is located inside of the center sphere conceal- 
ment at (—0.18m, 0.0, 0.0), the excited EM wave excied 
by S\ is named as First electric wave, E XXj i The sec- 
ond current point source S2 is located in free space at 
(0.73m, 0.0, 0.0) where is located the right side outside of 
the whole GL double layer cloaks. The EM wave excited 
by the S2 is named as Second electric wave.E xx 2, 
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FIG. 1: (color online) At time step 48dt, front of 
Second EM wave, E XX] 2 inside of the outer layer of GLWF 
cloak R2 < r < R3 propagates no faster than light speed, 
it is faster than the wave speed in figure 2 and in figure 3. 
The wave front of the First electric wave , E xx ,i, propagates 
inside the inner layer, Ri < r < Ri. 
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C. Comparison Between EM Wave Propagation 
Through The GLWF, GLO,And GLI-PSO Double 
Layer Cloaks 

The GL modeling simulations of the EM wave ex- 
cited by above two point sources Si and 5*2 propaga- 
tion through the GLWF, GLO,and GLI-PSO double layer 
cloaks are presented in the Figures 1-3 at 48th time step 
, in Figures 4-6 at 74th time step, in Figures 7-9 at 128th 
time step, respectively. For comparison, we arrange Fig- 
ures 1-3, Figures 4-6 , and Figures 7-9 as four figure group 
rows. At time step 48dt, front of electric wave E xx> 2 in- 
side of the outer layer of GLWF, GLO outer layer cloak, 
and PS outer layer, i?2 < r < R3 are presented in Figures 
1 -3 respectively. In Figure 1, E xx ,2 inside of the outer 
layer of GLWF, R 2 < r < R 3 , propagates no faster than 
light speed, it is faster than the wave speed in figure 2 
and in figure 3. In Figure 2 electric wave E xx ^ inside 
of the outer layer of GL outer layer cloak R2 < r < R3 
propagates no faster than light speed. It is slower than 
E xx ,2 in figure 1 and faster than E xx ^, in figure 3. In Fig- 
ure 3, front of E xx ^ inside of the outer layer of PS outer 
layer cloak, R2 < r < R3 , propagates no faster than light 
speed. It is slower than E xx ,2 in figure 1 and figure 2. 
The wave front of the E XXi ±, propagates inside the inner 
layer, Ri < r < i?2- At time step 74dt, front of elec- 
tric wave E XX: 2 inside of the outer layer of GLWF, GLO 
outer layer cloak, and PS outer layer, R 2 < r < R3 are 
presented in Figures 4 -6 respectively. In Figure 4, front 
of E xx ^2 inside of GLWF outer layer cloak R2 < r < R3 
propagates no faster than light speed, it is slower than 
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FIG. 2: (color online) At time step 48dt, front of 
Second EM wave, E XXt 2 inside of the outer layer of GL outer 
layer cloak R2 < r < R3 propagates no faster than light speed. 
It is slower than E xx ,2 in figure 1 and faster than E XXt 2 in fig- 
ure 3 The wave front of the E X x,x, propagates inside the inner 
layer, R\ < r < i?2- 



the wave speed in figure 5 and in figure 6. In Figure 5, 
front of E xx 2 inside of GL outer layer cloak R 2 < r < R 3 
propagates little faster than light speed. It is faster than 
E xx ,2 in GLWF in figure 4 and slower than E xx p in PS 
in figure 6. In Figure 6, front of E XXy 2 inside of PS outer 
layer cloak i?2 < r < R3 propagates faster than light 
speed. It is faster than E XX: 2 in GLWF in figure 4 and 
GLO in figure 5. At time step 128dt, front of electric wave 
E xx . 2 inside of the outer layer of GLWF, GLO outer layer 
cloak, and PS outer layer, R 2 < r < R 3 are presented in 
Figures 7 -9 respectively. In Figure 7, front of E xx ^ in- 
side of GLWF outer layer cloak R 2 < r < R 3 propagates 
slower than light speed, it is slower than the wave speed 
in figure 8 and in figure 9. In Figure 8, front of E XXi 2 
inside of GL outer layer cloak R2 < r < R3 propagates 
little faster than light speed. It is faster than E xx ^ in 
GLWF in figure 7 and slower than E xx $ in PS in figure 9. 
In Figure 9, front of E xx ^ inside of PS outer layer cloak 
i?2 < r < R3 propagates more faster than light speed. 
It is faster than E xx ^ in GLWF in figure 7 and GLO in 
figure 8. The wave front of the E xx ^i, propagates inside 
the inner layer, Ri < r < i?2- 



IV. THEORY OF RECIPROCAL LAW OF THE 
EM WAVE FIELD THROUGH THE CLOAKS 

A. Theory Of The EM Wave Field Through The 
GL Double Layer Cloaks 

We propose the theoretical analysis of the interaction 
between the EM wave and GL cloaks in this section. 
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FIG. 3: (color online) At time step 48dt, front of 
Second EM wave, E XXi 2 inside of the outer layer of PS outer 
layer cloak R2 < r < Rs propagates no faster than light speed. 
It is slower than E XXi 2 in figure 1 and figure 2. The wave front 
of the E XXt i, propagates inside the inner layer, R\ < r < R2- 
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FIG. 4: (color online) At time step 74dt, front of E XXi 2 inside 
of GLWF outer layer cloak R2 < r < R3 propagates no faster 
than light speed, it is slower than the wave speed in figure 5 
and in figure 6. The wave front of the First electric wave , 
E XXt i, propagates inside the inner layer, Ri < r < R2. 



FIG. 6: (color online)At time step 74dt, front of E xx> % inside 
of PS outer layer cloak R2 < r < R3 propagates faster than 
light speed. It is faster than E XXt 2 in GLWF in figure 4 and 
GLO in figure 5. The wave front of the E Km ,i, propagates 
inside the inner layer, Ri < r < Ri. 




FIG. 5: (color online) At time step 74dt, front of E xx ,2 inside 
of GL outer layer cloak R2 < r < R3 propagates a little faster 
than light speed. It is faster than E xx ,2 in GLWF in figure 4 
and slower than E XXt 2 in PS in figure 6 The wave front of the 
E xx ,i, propagates inside the inner layer, Ri < r < ife- 



FIG. 7: (color online) At time step 128dt, front of E XX} 2 inside 
of GLWF outer layer cloak R2 < r < R3 propagates slower 
than light speed, it is slower than the wave speed in figure 8 
and in figure 9. The wave front of the First electric wave , 
E XXi i, propagates inside the inner layer, Ri < r < .Re- 



statement 1: Let domain floL in (3) and the meta- 
material Dql in (4) be GL double layer cloak, and 
e=£b,M = Mb be basic permittivity and permeability, 
respectively, inside of the central sphere concealment 
|r'| < i?i and outside of the GL cloak \r'\ > i? 3 , we have 
the following statements: (1) provide the local source 
is located inside of the concealment of GL double layer 



cloak, |r s | < Ri, the excited EM wave field inside of the 
concealment never be disturbed by the cloak; (2) provide 
the local source is located inside of concealment or inside 
of the inner layer of the GL double layer cloak, \r s \ < R2, 
the EM wave field is vanished outside of the inner layer 
of GL cloak and is always propagating and going to the 
boundary r = R2 and before r — i?2- (3) provide the 
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FIG. 8: (color online) At time step 128dt, front of E XXy 2 
inside of GL outer layer cloak R2 <r < Rz propagates faster 
than light speed. It is faster than E XXi 2 in GLWF in figure 7 
and slower than E XXi 2 in PS in figure 9 The wave front of the 
E XXl i, propagates inside the inner layer, Ri < r < Ri- 
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FIG. 10: (color online) A double cloth anti detection is 
around the fly model, The figure 5 is from the figure 11 in 
paper [8] in 2001 
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FIG. 11: (color online) A double cloth anti detection is 
around the bar model, The figure 6 is from the figure 2 in 
paper [8] in 2001 



FIG. 9: (color online) At time step 128di, front of E xx ,2 
inside of PS outer layer cloak R2 < r < Rz propagates more 
faster than light speed. It is faster than E XXi 2 in GLWF in 
figure 7 and GLO in figure 8. The wave front of the E XXt i, 
propagates inside the inner layer, Ri < r < -R2. 



source is located outside of the GL double layer cloak, 
\f s \ > i?3, the excited EM wave field propagation outside 
of the double layer cloak as same as in free space and 
never be disturbed by the double layer cloak; (4) provide 
the local source is located outside of double layer cloak or 
located inside of the outer layer of GL cloak, \r s \ > R 2 , 
the excited EM wave field never propagate into the inner 
layer of GL cloak and the concealment. 



B. Theory Of Reciprocal Law Of The EM Wave 
Field Through The Cloaks 

Statement 2: (1) In the domain consist of free space, 
single layer cloak and its cloaked concealment with nor- 
mal material, the two sources reciprocal law is damaged. 
(2) In the domain consist of free space, single layer cloak 
and its cloaked concealment with some special double 
negative refractive index metamaterial, the two sources 
reciprocal law is recovered, but the cloak invisibility func- 
tion is lose. (3) In the domain consist of free space, GL 
double layer cloak and its cloaked concealment with nor- 
mal material, the two sources reciprocal law is satisfied, 
and the cloak invisibility function is complete and suffi- 
cient in wide frequency band. 
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C. There Exists No Maxwell EM Wavefield Can 
Be Excited By Nonzero Local Sources Inside Of The 
The Single Layer Cloaked Concealment With 
Normal Materials 



Statement 3: Suppose that a 3D anisotropic inhomo- 
geneous single layer cloak domain separates the whole 3D 
space into three sub domains, one is the single layer cloak 
domain fl c ik with the cloak material; the second one is 
the cloaked concealment domain fl CO ni with normal EM 
materials; other one is the free space outside of the cloak. 
If the Maxwell EM wavefield excited by a point source 
or local sources outside of the concealment f2 con / is van- 
ished inside of the concealment £l CO ni , then there excists 
no Maxwell EM wave field can be excited by the local 
sources inside the cloaked concealment £l CO ni with nor- 
mal materials. 

The statement 2 is proved by the GL method in au- 
thor's paper [12]. 



V. HISTORY AND DISCUSSIONS 
A. History 

A double layer cloth phenomenon to prevent the 
GILD inversion [6] [8] detection has been observed 
in paper [9] in 2001 which is published in SEG online 



B. Advantages Of The GL Method 

The GL EM modeling is fully different from FEM and 
FD and Born approximation methods and overcome their 
difficulties. There is no big matrix equation to solve in 
GL method. Moreover, it does not need artificial bound- 
ary and absorption condition to truncate the infinite do- 
main. Born Approximation is a conventional method in 
the quantum mechanics and solid physics However , it 
is one iteration only in whole domain which is not accu- 
rate for high frequency and for high contrast materials. 
The GL method divides the domain as a set of small sub 
domains or sub lattices. The Global field is updated by 
the local field from the interaction between the global 
field and local subdomain materials successively. Once 
all subdomain materials are scattered, the GL field solu- 
tion is obtained which is much more accurate than the 
Born approximation. GL method is suitable for all fre- 
quency and high contrast materials. When the size of 
the sub domain is going to zero, the GL method is con- 
vergent and has 0(h 2 ) if the trapezoidal integral formula 
is used, moreover, is has super convergence 0(ft. 4 ) if the 
Gaussian integral formula is used[10]. Chen et al pro- 
posed an analytical method for analysis of the PS cloak 
[7] . The GL method has double capabilities of the theo- 
retical analysis and numerical simulations that has been 
shown in this paper. 



|http: / / segdl.org/journals /doc/ S EGLIB-home /dci / searchDCI.jsp 
The double layer cloth to cloak fly from the exterior 
wave GILD detection is obvious around the fly which 
is shown in figure 10; the double cloth around the 
bar is shown in figure 11. We developed a novel 
and effective Global and Local field (GL) modeling 
and inversion[2] [3] [5] to study the meta materials, 
periodic photonic crystals and condense physics etc. 
wide physical sciences. 3D GL EM modeling and 
inversion [3] [5] and computational mirage have been 
presented in PIERS 2005 and published in proceeding of 
PIERS 2005 in Hangzhou, which can be downloaded from 
http://piers.mit.edu/piersproceedings/piers2k5Proc.php, 
please see the references of [2J. We developed 3D FEM 
for the elastic mechanics first in China in 1972 [10] and 
discovered the superconvergence of the 3D cubic curve 
isoparameter element first in the world [11]. The 3D 
isoparameter element can be used for making arbitrary 
curve cloak [10]. We deeply to know the merits and 
drawbacks of FEM. The GL method overcomes the 
drawbacks of FEM and FD methods. The history of 
development of FEM and GILD and GL method has 
been described in [11] and reference of [2]. The 3D 
and 2D GL parallel software is made and patented by 
GLGEO. The GL modeling and its inversion [3] [5] and 
GL EM quantum field modeling are suitable to solve 
quantization scattering problem of the electromagnetic 
field in the dispersive and loss metamaterials, cloaks and 
more wide anisotropic materials. 



VI. CONCLUSIONS 



The simulations of the EM wave propagation through 
the GLWF. GLO, and GLI-PSO double layer cloaks and 
comparison between them show that the GLWF and 
GLO double layer cloak overcomes the following diffi- 
culties of the single layer PS cloak. (1) The PS cloak 
damaged the EM environment of its concealment, such 
that there exists no EM wavefield can be excited inside 
the concealment of the PS cloak, the concealment of the 
PS cloak is blind. Our GL double layer cloak recovered 
the normal EM environment in its concealment, such that 
the EM wave field can be excited inside the concealment 
of the GL double layer cloak. (2) The PS cloak is very 
strong dispersive and strong degenerative cloak material. 
The PS cloak has invisibility only in very narrow fre- 
quency band. There is exceeding light speed physical 
violation in PS cloak. Our GLWF double layer cloak 
corrects the violation. (3) The two sources reciprocal 
law is very important principle in the electromagnetic 
theory and application. The reciprocal law is satisfied 
in our GLWF and GLO double cloak media. However, 
the PS cloak damaged the reciprocal law. The follow- 
ing physical statements are described: (1) In the domain 
consist of free space, single layer PS cloak and its cloaked 
concealment with normal material, the two sources recip- 
rocal law is damaged. (2) In the domain consist of free 
space, single layer cloak and its cloaked concealment with 
some special double negative refractive index metamate- 
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rial, the two sources reciprocal law is recovered, but the 
cloak invisibility function is lose. (3) In the domain con- 
sist of free space, GLWF and GLO double layer cloak 
and its cloaked concealment with any material, the two 
sources reciprocal law is satisfied, and the cloak invisibil- 
ity function is complete, sufficient, moreover the GLWF 
has all advantages of the GL double layer cloak in broad 
frequency band. The GLWF and GLO double layer cloak 
materials and the 3D and 2D GL parallel algorithms and 
software are made by authors in GL Geophysical Lab- 
oratory and are patented by GLGEO and all rights are 
reserved in GLGEO. 

The GL method is an effective physical simulation 
method. It has double abilities of the theoretical analysis 



and numerical simulations to study the cloak metamate- 
rials and wide material and Field scattering in physical 
sciences. 
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